A=[1,0,0,0;1.1,1,0,0;1.21,1.1,1,0;1.331,1.21,1.1,1];
b=[400;440;484;532.4];
x0=[100;100;100;100];
lb=0;
options=optimset('Display','off');
[x,f]=fmincon(@f7,x0,A,b,[],[],lb,[],[],options)
-f